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^ Abstract: Experimental data from the National Air Surveillance Network of Japan from 

t;;}^ 1974 to 1996 and from independent measurements performed simultaneously in the regions 

2 of Ljubljana (Slovenia), Odessa (Ukraine) and the Ukrainian "Academician Yernadsky" 

^ Antarctic station (64°15W; 65°15S), where the air elemental composition was determined 

by the standard method of atmospheric particulate matter (PM) collection on nucleopore 
>< m,e. and ..b.e,.e„t neu.cn acivaUon analy.i. were analyzed. Con^para.ive anal,., of 

different pairs of atmospheric PM element concentration data sets, measured in different re- 
gions of the Earth, revealed a stable linear (on a logarithmic scale) correlation, showing a 
power law increase of every atmospheric PM element mass and simultaneously the cause of 
this increase - fractal nature of atmospheric PM genesis. Within the framework of multi- 
fractal geometry we show that the mass (volume) distribution of atmospheric PM elemental 
components is a log normal distribution, which on a logarithmic scale with respect to the 
random variable (elemental component mass) is identical to normal distribution. This means 
that the parameters of two-dimensional normal distribution with respect to corresponding at- 
mospheric PM-multifractal elemental components measured in different regions, are a priory 
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connected by equations of direct and inverse linear regression, and the experimental manifes- 
tation of this fact is the linear correlation between the concentrations of the same elemental 
components in different sets of experimental atmospheric PM data. 

Keywords: Atmospheric aerosols; Multifractal; Neutron activation analysis; South Pole; 
Ukrainian Antarctic station 



1. Introduction 

Analysis of the concentrations of elements characteristics of the terrestrial crust, anthropogenic emis- 
sions and marine elements used in monitoring of the levels of atmospheric aerosol contamination in- 
dicates that these levels at the two extremes of: (i) the Antarctic (South Pole [1]), and (ii) the global 
constituent of atmospheric contamination measured at continental background stations [2] display sim- 
ilar patterns. A distinction between them is, however, evident and consists in that the mean element 
concentrations in the atmosphere over continental background stations (Ccb) located in different re- 
gions of the Earth exceed the corresponding concentrations at the South Pole (Csp) by some 20 — 10^ 
times. 

Comparing the concentrations of a given element / in atmospheric aerosol from samples {Csp,i} and 
{CcB,i}, it became evident that the dependence of the mean concentrations of any particular element 
from the sample {Csp,i} or {CcB,i} on a logarithmic scale is described by a linear one, with good 
precision for any of the elements from a given pair of sampling station: 

^^C^.^ = acB-sp + bcB~sp^^Cgp, (1) 

bcB-sp ~ 1, (2) 

where acs-sp and hcs-sp are the intercept and slope (regression coefficient) of the regression line, 
respectively. 

This unique and rather unexpected result was first established by Pushkin and Mikhailov [2]. It is 
noteworthy, according to [2], that the reason for the large enrichment of atmospheric aerosol with those 
elements which are exceptions to the linear dependence, is related mostly either to the anthropogenic 
contribution produced by extensive technological activity, e.g. the toxic elements (Sb, Pb, Zn, Cd, As, 
Hg), or to the nearby sea or ocean or sea as a powerful source of marine aerosol components (Na, I, Br, 
Se, S, Hg). 

However, our numerous experimental data specify persistently that the Pushkin-Mikhailov depen- 
dence (1) in actual fact is the particular case (at hi2 = 1) of more general of linear regression equation 

In (Cii/pi) = ai2 + hi2 In {€21/ Pi) (3) 

where Ci and pi are the concentration and specific density of i-th isotope component in atmospheric PM 
measured in different regions (the indexes 1 and 2) of the Earth. 
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It is obvious, if we will be able to prove that the linear relation (3) in element concentrations between 
the above mentioned samples reflects the more general fundamental dependence, it can be used for the 
theoretical and experimental comparison of atmospheric PM independently of the given region of the 
Earth. Moreover, the linear relation (3) can become a good indicator of the elements defining the level 
of atmospheric anthropogenic pollution, and thereby to become the basis of method for determining a 
pure air standard or, to put it otherwise, the standard of the natural level of atmospheric pollution of 
different suburban zones. This is also indicated by the power law character of Eq. (3), reflecting the fact 
that the total genesis of non-anthropological (i.e natural) atmospheric aerosols does not depend upon the 
geography of their origin and is of a fractal nature. 

In our opinion, this does not contradict the existing concepts of microphysics of aerosol creation and 
evolution [3], if we consider the fractal structure of secondary {Dp > 1/im) aerosols as structures formed 
on the prime inoculating centres, (Dp < 1/im). Such a division of aerosols into two classes - primary and 
secondary [3] - is very important since it plays the key role for understanding of the fractal mechanism 
of secondary aerosol formation, which show scaling structure with well-defined typical scales during 
aggregation on inoculating centres (primary aerosols) [4]. 

The objective of this work was twofold: (i) to prove reliably the linear validity of Eq. (3) through 
independent measurements with good statistics performed at different latitudes and (ii) to substantiate 
theoretically and to expose the fractal mechanism of interrelation between the genesis, size and com- 
position of atmospheric PM measured in different regions of the Earth, in particular in the vicinity of 
Odessa (Ukraine), Ljubljana (Slovenia) and the Ukrainian Antarctic station "Academician Vernadsky" 
(64°15W; 65°15S). 

2. The linear regression equation and experimental data of National Air Surveillance Network 
of Japan 

In this study, experimental data [5] from the National Air Surveillance Network (NASN) of Japan 
for selected crustal elements (Al, Ca, Fe, Mn, Sc and Ti), anthropogenic elements (As, Cu, Cr, Ni, 
Pb, V and Zn) and a marine element (Na) in atmospheric particulate matter obtained in Japan for 23 
years from 1974 to 1996 were evaluated. NASN operated 16 sampling stations (Nopporo, Sapporo, 
Nonotake, Sendai, Niigata, Tokyo, Kawasaki, Nagoya, Kyoto-Hachiman, Osaka, Amagasaki, Kurashiki, 
Matsue, Ube, Chikugo-Ogori and Ohmuta) in Japan, at which atmospheric PM were regularly collected 
every month by a low volume air sampler and analyzed by neutron activation analysis (NAA) and X-ray 
fluorescence (XRF). During the evaluation, the annual average concentration of each element based on 
12 monthly averaged data between April (beginning of financial year in Japan) and March was taken 
from NASN data reports. The long-term (23 years) average concentrations were determined from the 
annual average concentration of each element. 

Analysis of the NASN data [5] shows that the highest average concentrations were observed in 
Kawasaki (Fe, Ti, Mn, Cu, Ni and V), Osaka (Na, Cr, Pb and Zn), Ohmuta (Ca) and Niigata (As), 
respectively. These cities are either industrial or large cities of Japan. Conversely, the lowest average 
concentrations were noticed in Nonotake (Al, Ca, Fe, Ti, Cu, Cr, Ni, V and Zn) and Nopporo (Mn, As 
and Pb), as expected [5]. On the basis of these results, Nonotake and Nopporo were selected as the 
baseline-remote area in Japan. 
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A simple model of linear regression, in which the evaluations were made by the least squares method 
[6,7], was used to build the linear dependence described by Eq. (3). The results of NASN data presented 
on a logarithmic scale relative to the data of Nonotake {C*AronotaA;e} (^^^ Figure 1) show with high con- 
fidence the adequacy of experimental and theoretical dependence of Eq. (3) type. As can be seen from 
Figure 1, the Nonotake station was chosen as the baseline, where the lowest concentrations of crustal 
and anthropogenic elements and small variations in time (23 years) were observed [5]. The NASN data 
are presented in the form of "city - city" concentration dependences. 

Figure 1. Relation between annual concentrations of elements in atmospheric particulate 
matter over Japan and the same data obtained in the region of Nonotake (data from NASN 
of Japan [5]). In some cases (see text) the concentration of anthropogenic element Cr was 
excluded from the data. The underlined cities are large industrial centres in the Japan [5]. 
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Analysis of concentration data for Japanese city atmospheric PM unambiguously shows that the i-th 
element mass in atmospheric PM grows by power low, proving the assumption [4,8-10] about the fractal 
nature of atmospheric PM genesis. 

3. The linear regression equation and the composition of atmospheric aerosols in different re- 
gions of the Earth 

It is evident that in order to generalize the results from NASN data processing more widely, the valid- 
ity of Eq. (3) should be checked on the basis of atmospheric aerosol studies in performed independent 
experiments at different latitudes. For this reason such studies were performed in the regions of Odessa 
(Ukraine), Ljubljana (Slovenia) and the Ukrainian Academecian Vemadsky Antarctic station (64°15W; 
65°15S). The determination of the element composition of the atmospheric air in these experiments was 
performed by the traditional method based on collection of atmospheric aerosol particles on nucleopore 
filters with subsequent use of fco-instrumental neutron activation analysis. Regression analysis was used 
for processing of the experimental data. 
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Figure 2. Schematic representation of the flow-chamber of the PMIO type device for air 
sampling. 




ISmin 



3.1. Experimental 

3.1.1. Collection of atmospheric aerosol particles on nucleopore filters 

For collection of atmospheric aerosol particles on filters, a device of the PMIO type was used with 
the Gent Stacked Filter Unit (SFU) interface [11,12]. The main part of this device is a flow-chamber 
containing an impactor, the throughput capacity of which is equivalent to the action of a filter with an 
aerodynamic pore diameter of 10 /im and a 50% aerosol particle collection efficiency based on mass, 
and an SFU interface designed by the Norwegian Institute for Air Research (NILU) and containing two 
filters (Nucleopore) each 47 mm in diameter, the first filter with a pore diameter of 8 /xm and the second 
filter with 0.4 fim pore diameter. It was experimentally found [12] that such geometry (Figure 2) results 
in an aerosol collection efficiency of the first filter of approximately 50%, whereas for the second filter 
this value was close to 100% [12]. More detailed results of thorough testing of a similar device can be 
found in [12]. 

3.1.2. ko-instrumental neutron activation analysis 

Airborne particulate matter (APM) loaded filters were pelleted with a manual press to a pellet of 5 
mm diameter and each packed in a polyethylene ampoule, together with an Al-0.1%Au IRMM-530 disk 
6 mm in diameter and 0.2 mm thickness and irradiated for short irradiations (2-5 min) in the pneumatic 
tube (PT) of the 250 kW TRIGA Mark II reactor of the J. Stefan Institute at a thermal neutron flux of 
3.5 ■ 10^^cm~^s~^, and for longer irradiations in the carousel facility (CF) at a thermal neutron flux of 1 . 1 ■ 
10^^cm~^s~^ (irradiation time for each sample about 18-20 h). After irradiation, the sample and standard 
were transferred to clean 5 mL polyethylene mini scintillation vial for gamma ray measurement. To 
determine the ratio of the thermal to epithermal neutron flux (/) and the parameter a, which characterizes 
the degree of deviation of the epithermal neutron flux from the 1 / E-\aw, the cadmium ratio method for 
multi monitor was used [13]. It was found that / = 32.9 and a = —0.026 in the case of the PT channel, 
and / = 28.7 and a = —0.015 for the CF channel. These values were used in calculation of the 
concentrations of short- and long-lived nuclides. 
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7-activity of irradiated samples were measured on two HPGe-detectors (ORTEC, USA) of 20 and 
40% measurement efficiency [13]. Experimental data obtained on these detectors were fed into and pro- 
cessed on EG&G ORTEC Spectrum Master and Canberra SlOO high-velocity multichannel analyzers, 
respectively. To calculate net peak areas, HYPERMET-PC V5.0 software was used [14], whereas for 
evaluation of elemental concentrations in atmospheric aerosol particles, KAYZERO/SOLCOI® soft- 
ware was used [15]. More details of the fco-instrumental neutron- activation analysis applied could be 
found in [13]. 

3.2. Comparative analysis of atmospheric PM composition in different regions of the Earth 

The results of presenting the atmospheric PM concentration values of Ljubljana {C\^juhijana\ ^'^d 
Odessa {CQdessa} ^ logarithmic scale relative to the similar data from Academician Vemadsky station 
{^Ant station} dcmoustratc with high reliability that the correlation coefficient r is approximately equal to 
unity both for the direct and reverse regression lines "Odessa - Antarctic station", "Ljubljana - Antarctic 
station" (Figure 3): 

r = [6i2&2i]'^' ~ 1, (4) 

where and 621 are the slopes of direct and reverse regression lines (Eq. (1)) for corresponding pairs. 

Figures 4 and 5 show the regression hnes for daily normalized average concentrations of crustal, an- 
thropogenic and marine elements (Table 1), measured on March 2002 in the regions of Odessa (Ukraine), 
the Ukrainian Antarctic station (°15W; 65°15S) and Ljubljana (Slovenia) relative to the same data ob- 
tained in Nonotake [5] (Figure 4) and the South Pole [1] (Figure 5). 

Table 1. Chemical element composition in atmospheric PM. Symbols Csp, Ccb, C Nonotake 
denote concentrations at the South Pole, continental background stations and Nonotake-city, 
respectively. Measured concentrations in the vicinity of the Ukrainian Antarctic Station, 
Odessa and Ljubljana are denoted by C Antarctica, Codessa and Cijubijana 



Element 


Csp 


Ccb 


C Nonotake 


C Antarctica 


C Odessa 


C Ljubljana 


{ng/m?) 


S 


49 












Si 














CI 


2.4 


90 








528.2 


Al 


0.82 


1.2 X 10^ 








1152 


Ca 


0.49 




185.8 


93 


2230 


1630 


Fe 


0.62 


1 X 10 2 


157.6 


53.4 


1201 


1532 


Mg 


0.72 










1264 


K 


0.68 






24.4 


502.8 


918 


Na 


3.3 


1.4 X 10 2 


538.8 


361.95 


393.8 


1046 


Pb 




10 


23.2 








Zn 


3.3 X 10-2 


10 


38.3 


4.48 


62.6 


127.2 
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Ti 


0.1 




15.4 








F 














Br 


2.6 


4 




1.34 


14.81 


9.62 


Cu 


2 X 10 


3 


8.20 


11050 


4240 




Mn 


1.2 X 10 


3 


6.61 






34.7 


Ni 




1 


1.48 








Ba 


1.6 X 10 










6.91 


V 


1.3 X 10 


1 


2.44 






17.08 


I 


0.74 










36.24 


Cr 


4 X 10 


0.8 


1.14 


3.11 


45.2 




Sr 


5.2 X 10 












As 


3.1 X 10 


1 


2.74 




1.83 


4.18 


Rb 


2 X 10 










0.92 


Sb 


8 X 10 


0.5 




0.07 


1.97 


4.62 


Cd 


< 1.5 X 10 


0.4 










Mo 


49 






1.46 


4.08 




Se 


< 0.8 


0.3 




0.02 


0.47 


1.21 


Ce 


4 X 10 








4.22 


2.72 




0.17 


0.3 




0.36 


2.75 




W 


1.5 X 10 






0.24 






La 


4.5 X 10 ~^ 








1.79 


1.36 


Ga 


< 1 X 10 












Co 


5 X 10 


0.1 




0.02 


0.89 


0.29 




< 4 X 10 






2.74 


1.42 




Cs 


1 X 10 








0.09 




So 


1.6 X 10 


5 X 10 


0.04 


3 36 X 10 


0.21 


0.298 


Th 


1.4 X 10 








0.29 


0.060 


u 










0.07 




Sm 


9 X 10 






2 68 X 10 


0.24 


0.198 


In 


5 X 10 












Ta 


7 X 10 












Hf 


6 X 10-5 












Yb 


<0.05 












Eu 


2 X 10-5 












Au 


4 X 10 -5 






1.07 X 10-5 


0.006 




Lu 


6.7 X 10-6 













The choice of the atmospheric aerosol concentration values of the South Pole as a baseline, relative 
to which a dependence of the type given by Eq. (1) was analyzed, was made for three reasons: (i) these 
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Figure 3. Lines of direct regression "Odessa - Antarctic station" (a), "Ljubljana - Antarctic 
Station" (c) and line of inverse regression "Antarctic Station - Odessa" (b), "Antarctic Station 
- Ljubljana" (d). 
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Figure 4. Relationship between atmospheric PM elemental concentrations measured in the 
regions Odessa (Ukraine), Ljubljana (Slovenia), Vemadsky station (64°15W; 65°15S), South 
Pole [1] and the same data measured in the region of Nonotake [5]. 
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Figure 5. Relationship between atmospheric PM elemental concentrations measured in 
the regions Odessa (Ukraine), Ljubljana (Slovenia), Vernadsky station (64°15W; 65°15S), 
Nonotake [5] and the same data measured in the region of the South Pole [1]. 




South Pole, g/crrf 
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data were obtained by the same technique and method as in section 3.1, and simultaneously expand 
the geographical comparison, (ii) the element spectrum that characterizes the atmosphere of South Pole 
covers a wide range of elements (see Table 1); and (iii) the South Pole has the purest atmosphere on 
Earth, making it a convenient basis for comparative analysis. 

We present also the monthly normalized average concentrations of atmospheric PM measured over 
the period 2006-2007 in the region of the Ukrainian Antarctic station "Academician Vernadsky" (Figures 
6, 7). 

Comparative analysis of experimental sets of normalized concentrations of atmospheric aerosol ele- 
ments measured in our experiments and independent experiments of the Japanese National Air Surveil- 
lance Network (NASN) shows a stable linear (on a logarithmic scale) dependence on different time scales 
(from average daily to annual). That points to a power law increase of every atmospheric PM element 
mass (volume) and simultaneously to the cause of this increase - the fractal nature of atmospheric PM 
genesis. 

In other words, stable fulfillment of the equality (4) not only for the experimental data shown in 
Figures 3-7, but also for any pairs of the NASN data [5] (see Figure 2) unambiguously indicates that, on 
the one hand, the model of linear regression is satisfactory and, on the other hand, any of the analyzed 
samples {m, }, which describe the sequence of i-th element partial concentrations in an aerosol, must 
to obey the Gauss distribution with respect to the random quantity \npi. Proof of these assertions for 
multifractal objects is presented below. 

4. The spectrum of multifractal dimensions and log normal mass distribution of secondary aerosol 
elements 

A detailed analysis of Figs. 1, 3-5, where the linear regressions for different pairs of experimental 
samples of element concentrations in atmospheric PM measured at various latitudes are shown, allows 
us to draw a definite conclusion about the multifractal nature of PM. The basis of such a conclusion is 
the reliably observed linear dependence of Eq. (3) type between the normalized concentrations Ci of the 
same element i in atmospheric PM in different regions of the Earth. 

Thus, it is necessary to consider the atmospheric PM, which is the multicomponent (with respect 
to elements) system, as a nonhomogeneous fractal object, i.e., as a multifractal. At the same time, 
the spectrum of fractal dimensions /(a) and not a single dimension ao (which is equal to Dq for a 
homogeneous fractal) is necessary for complete description of a nonhomogeneous fractal object. We 
will show below that the spectrum of fractal dimensions f{a) of multifractal predetermines the log 
normal type of statistics or, in other words, the log normal type of mass distribution of multifractal i-th 
component. This is very important, because the representation of mass distribution of atmospheric PM 
as the log normal distribution is predicted within the framework of the self-preserving theory [20,21] and 
is confirmed by numerous experiments at the same time [3]. 

To explain the main idea of derivation we give the basic notions and definitions of the theory of 
multifractals. Let us consider a fractal object which occupies some bounded region £ of size L in 
Euclidian space of dimension d. At some stage of its construction let it represents the set of >> 1 
points distributed somehow in this region. We divide this region £ into cubic cells of side e « L and 
volume e'^. We will take into consideration only the occupied cells, where at least one point is. Let i 
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Figure 6. The regression lines for the normalized monthly average concentrations of atmo- 
spheric PM measured in the region of the Ukrainian Antarctic station "Academician Vemad- 
sky" measured in August-December 2006 with respect to October 2006. 




October, g/cm ^ 
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Figure 7. The regression lines for the normalized monthly average concentrations of atmo- 
spheric PM measured in the region of the Ukrainian Antarctic station "Academician Vemad- 
sky" in January-March 2007 with respect to October 2006. 




October, g/cm ' 
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be the number of occupied cell i = 1, 2, N{e), where N{e) is the total number of occupied cells, 

which depends on the cell size e . Then in the case of a regular (homogeneous) fractal, according to the 
definition of fractal dimensions D, the total number of occupied cells N(e) at quite small e looks like 

Nie) ^ el'' = Lf . (5) 

where e is the cell size in L units, is the size of fractal object in e units. 

When a fractal is nonhomogeneous a situation becomes more complex, because, as was noted above, 
a multifractal is characterized by the spectrum of fractal dimensions f{a), i.e. by the set of probabilities 
Pi, which show the fractional population of cells e by which an initial set is covered. The less the cell 
size is, the less its population. For self-similarity sets the dependence pi on the cell size e is the power 
function 



M^) = ^ = ^r"' (6) 



1 

N~(e) 

where is a certain exponent which, generally speaking, is different for the different cells i. It is 
obvious, that for a regular (homogeneous) fractal all the indexes ctj in Eq. (6) are identical and equal to 
the fractal dimension D. 

We now pass on to probability distribution of the different values Oj. Let n{a)da is the probability 
what tti is in the interval [a, a + da]. In other words, n(a)da is the relative number of the cells i, which 
have the same measure Pi as in the interval [a, a + da]. According to (5), this number is proportional 
to the total number of cells N{e) ~ ej^^ for a monofractal, since all of «i are identical and equal to the 
fractal dimension D. 

However, this is not true for a multifractal. The different values of occur with probability char- 
acterized by the different (depending on a) values of the exponent f{a), and this probability inherently 
corresponds to a spectrum of fractal dimensions of the homogeneous subsets £q of initial set £: 

n(a)^£Z^("^ (7) 

Thus, from here a term "multifractal" becomes clear. It can be understood as a hierarchical joining of 
the different but homogeneous fractal subsets £a of initial set £, and each of these subsets has the own 
value of fractal dimension f{a). 

Now we show how the function f{a) predetermines the log normal kind of mass (volume) distribu- 
tion of the multifractal i-th component. To ease further description we represent expression (8) in the 
following equivalent form: 

n{a) ^ exp [f{a) In L^] . (8) 

It is not hard to show [19] that the single-mode function f(a) can be approximated by a parabola near 
its maximum at the value oq. 

f{a)^Do-v{a-aoy, (9) 
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where the curvature of parabola 



V 



r(«o) 



(10) 



2 2 [2(ao - ^o) + D'^^,] 
is determined by the second derivative of function /(a) at a point Uq. Due to a convexity of the function 
/(a) it is obvious, that the magnitude in square brackets must be always positive. The fact, that the last 
summand D'^^q in these brackets is numerically small and it can be neglected, will be grounded below. 
At the large the distribution n(a) (8) with an allowance for (9) takes on the form 



n{a) ~ exp 



{a — ao) lnZ/£ 



4 (ao - ^o) 

Then, taking into account (5), we obtain from (11) the distribution function of random variable p, 



(11) 



n{pi) InNoo ■ exp 



1 



41n(l/poiVBo) 

which with consideration of normalization takes on the final form 



(In Pi - InpoY 



(12) 



1 



27r(T 



exp 



-a o 

^ 2 



exp 



2^2 



{\npi-^f 



(13) 



where 



H = Inpo, cr^ = 2 In 



PoNdo 



(14) 



This is the so-called log normal (relative) mass Pi distribution. It is possible to present the first 
moments of such a kind of distribution for random variable pi in the following form: 



{p) = exp ( /i + -a' 



1/pIn' 



Do 



^2ao-3Do 



(15) 



var(p) = exp (2/i) ■ [exp (4^2) - exp (Sa^)] = Lf^^o-3Do) (^j^2{2ao-Do) _ _ 



(16) 



At the same time, it is easy to show that the distribution (13) for the random variable Inpi has the 
classical Gaussian form 



P(lnp) 



V2 



exp 



2a 



1 2 

[\np — (Inp)] 



(17) 



where the first moments of this distribution for the random variable Inp^ look like 



(Inp) =n + a^ = {ao- 2Do) In L, 



(18) 



var (Inp) = cr^ = 2 (ao — -Do) InL^. 



(19) 
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Thus, according to the known theorem of multidimensional normal distribution shape [22], a normal 
law of plane distribution for the two-dimensional random variable (pu, p2i) will be written down as 



P i\npii,\np2i) = ==exp 

/7r(Ticr2V 1 — r-' 



2(1 -r2) 



[u^ + f ^ — 2ruv) 



(20) 



where 



u = , v = , (21) 

r = cov(lnpij, lnp2i)/^i^2 is the correlation coefficient between Inpij and lnp2i- 

Then, by virtue of well-known linear correlation theorem [6,22] it is easy to show that hipu and lnp2i 
are connected by the linear correlation dependence, if the two-dimensional random variable (In pu , In p2i ) 
is normally distributed. This means that the parameters of two-dimensional normal distribution of the 
random values pu and p2i for the i-ih. component in one aerosol particle, which are measured in different 
regions of the Earth (the indexes 1 and 2), are connected by the equations of direct linear regression: 

Inpii - (Inpii) = r — [lnp2i - (lnp2i)] (22) 

0"2i 



and inverse linear regression 



lnp2i - (lnp2i) = r— [Inpij - (Inpij)] , (23) 
o"ii 



where i = 1, . . . ,Npis the component number. 

Taking into consideration that we measure experimentally the total concentration Cj of the z-th com- 
ponent in the unit volume of atmosphere, the partial concentration rrii of the i-th component in one 
aerosol particle measured in different regions of the Earth (the indexes 1 and 2) looks like 

mu = Cii/rii, m2i = C2i/n2, (24) 

where rii and n2 are the number of inoculating centers, whose role play the primary aerosols (Dp < 
1/im). 

Here it is necessary to make important digression concerning the choice of quantitative measure for 
description of fractal structures. According to Feder [18], determination of appropriate probabilities cor- 
responding to the chosen measure is the main difficulty. In other words, if choice of measure determines 
the search procedure of probabilities {pi}, which describe the increment of the chosen measure for given 
level of resolution e, then the probabilities themselves predetermine, in its turn, the proper method of 
their measurement. So, general strategy of quantitative description of fractal objects, in general case, 
should contain the following direct or reverse procedure: the choice of measure - the set of appropriate 
probabilities - the measuring method of these probabilities. 
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We choose the reverse procedure. So long as in the present work the averaged masses of elemental 
components of atmospheric PM-multifractal are measured, the geometrical probabilities, which can be 
constructed by experimental data for some fixed e, have the practically unambiguous form: 

Pi {e = const) = ^ — = =— — , (25) 

i i 

where pi is the specific gravity of the z-th component of secondary aerosol. 

Since the random nature of atmospheric PM formation is a priori determined by the random process 
of multicomponent diffusion-limited aggregation (DLA), we used the so-called harmonic measure [18] 
to describe quantitatively a stochastic surface inhomogeneity or, more precisely, to study an evolution of 
possible growth of cluster PM diameter. 

In practice, a harmonic measure is estimated in the following way. Because the perimeter of clusters, 
which form due to DLA, is proportional to their mass, the number of knots Np on the perimeter, i.e., 
the number of possible growing-points, is proportional to the number of cells in a cluster. Both 
these magnitudes, Np and A^, change according to the power law (5) depending on the cluster diameter 
L. From here it follows that all the knots Ap, which belong to the perimeter of such clusters, have a 
nonzero probability what a randomly wandering particle will turn out in them, i.e., they are the carriers 
of harmonic measure M(i{q, Sl)'- 



1=1 



oo, d < T{q) 



where Z(q, L^) is the generalized statistical sum in the interval — oo < g < oo, r(g) is the index of mass, 
at which a measure does not become zero or infinity at L — )■ oo(£:l — ^ 0). 

It is obvious, that in such a form the harmonic measure is described by the full index sequence 
r(g), which determines according to what power law the probabilities pi change depending on L. At 
the same time, the spectrum of fractal dimensions for the harmonic measure is calculated in the usual 
way, but using the "Brownian particles-probes" of fixed diameter e for study of possible growth of the 
cluster diameter L. From (26) it follows that in this case the generalized statistical sum Z(q, El) can be 
represented in the form 

Ziq,e,) = J^p^r^el'-^'\ (27) 

i=l 

As is known from numerical simulation of a harmonic measure, when the DLA cluster surface is 
probed by the large number of randomly wandering particles, the peaks of "high" asperities in such a 
fractal aggregate have greater possibilities than the peaks of "low" asperities. So, if possible growing- 
points on the perimeter of our aerosol multifractal to renumber by the index i = 1, . . . , Ap, the set of 
probabilities 

^ = {p^}i^ , (28) 
composed of the probabilities of Eq. (25) type will emulate the possible set of interaction cross-sections 
between Brownian particle and atmospheric PM-multifractal surface, which consists of the Ap groups 
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of identical atoms distributed on the surface. Each of these groups characterizes the i-th elemental 
component in the one atmospheric PM. 

A situation is intensified by the fact that by virtue of (17) each of the independent components obeys 
the Gauss distribution which, as is known [22], belong to the class of infinitely divisible distributions or, 
more specifically, to the class of so-called a-stable distributions. This means that although the Gauss dis- 
tribution has different parameters (the average fii = (Inpj) and variance af = var(lnpj) for each of com- 
ponents, the final distribution is the Gauss distribution too, but with the parameters fi = {Inp) = J2 l^i 
and CT^ = Yl '^i- From here it follows that the parameters of the two-dimensional normal distribution 
of all corresponding components in the plane {^1,^2} are connected by the equations of direct linear 
regression: 



Inpi - (Inpi) = r— [lnp2 - (lnp2)] (29) 
0-2 



and inverse linear regression: 



0"2 

\np2 - (lnp2) = r — [Inpi - (Inpi)] . (30) 

0"! 

So, validity of the assumption (28) for the secondary aerosol or, that is the same, validity of the 
application of harmonic measure for the quantitative description of aerosol stochastic fractal surface, 
will be proven if simultaneous consideration of (28) and equation of direct linear regression (29) will 
result in an equation identical to the equation of direct linear regression (3). 

Therefore, taking into account (28) we write down, for example, the equation of direct linear re- 
gression or, in other words, the condition of linear correlation between the samples of i-th component 
concentrations {In (Cu/pi)} and {In (C2i/pi)} in an atmospheric aerosol measured in different places 
(indexes 1 and 2): 

Np 

1 n ^^''^p^'^ 

\n{Cii/pi) = ai2 + bi2^n{C2i/pi), ai2 = — In ' ^ —7-, bi2 = r—. (31) 

Np ( Np X"^^ 02 

n {C2.IP. 

i=l 

It is obvious, that this equation completely coincides with the equation of linear regression (3), but 
is theoretically obtained on basis of the Gauss distribution of the random magnitude \npi and not in an 
empirical way. 

Physical interpretation of the intercept ai2 is evident from the expression (29), whereas meaning of 
the regression coefficient 612 becomes clear, according to (19) and (29), from the following expression: 



var (Inpi) 



bi2 = r { r = r 

var (In p2j 



1/2 



ln(Li/£i) 



1/2 



(32) 



_\n (12/62) _ 

where Li and L2 are the average sizes of separate atmospheric PM-multifractals typical for the atmo- 
sphere of investigated regions (indexes 1 and 2) of the Earth, ei and £2 are the cell sizes into which the 
corresponding atmospheric PM-multifractals are divided. 
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Below we give a computational procedure algorithm for identification of the generalized fractal di- 
mension Dq spectra and function /(a). It is obvious, that such a problem can be solved by the following 
redundant system of nonlinear equations of Eqs. (15), (16), (18) and (19) type: 

ln(p) = (2ao-3Do)ln(L/£), 

(33) 



var(p) = (L/£)^^^"''-^^o^ (L/e)^^^""-^"^ - 1 
(Inp) = /i + (T^ = (ao - 2Do) In (L/e), 

(34) 

var (Inp) =a^ = 2{ao- Dq) In (L/e), 

where e is the cubic cell fixed size, into which the bounded region £ of size L in Euclidian space of 
dimension d is divided. 

To solve the system of equations (33)-(34) with respect to the variables ao. Do and e, Sl it is necessary 
and sufficient to measure experimentally the i-th components of the concentration sample {Cj} in the 
unit volume of atmospheric air (see section 3) and size distribution of atmospheric PM for determination 
of the average size L. 

It will be recalled that from the physical standpoint so-called the box counting dimension Dq, the 
entropy dimension Di and the correlation dimension D2 are the most interesting in the spectrum of the 
generalized fractal dimensions Dg corresponding to different multifractal inhomogeneities. Within the 
framework of the notions and definitions of multifractal theory mentioned above we describe below the 
simple procedure for finding of spectrum of the generalized fractal dimensions Dg, taking into account 
the solution of the system of equations (33)-(34). 

From (27) it follows that in our case a multifractal is characterized by the nonlinear function r(g) of 
moments q 

In Z{q, El) 

r g = hm — . (35) 

As well as before we consider a fractal object, which occupies some bounded region £ of "running" 
size L (so that El 0) in Euclidian space of dimension d. Then spectrum of the generalized fractal 
dimensions Dg characterizing the multifractal statistical inhomogeneity (the distribution of points in the 
region £) is determined by the relation 

Dg = (36) 
q-1 

where (g — 1) is the numerical factor, which normalizes the function r(g) so that the equality Dg = dh 
fulfilled for a set of constant density in the rf-dimensional Euclidian space. 

Further, we are interested in the known in theory of multifractal connection between the mass index 
r(g) and the multifractal function /(a) by which the spectrum of generalized fractal dimensions Dg is 
determined 

= ^ = gr^I ■ ^(^) - f (^(^))] • (37) 



Version January 18, 2013 submitted to Atmosphere 



23 of 27 



It is obvious, that in our case, when the sample {pi} is experimentally determined and the cell size 
e « L is numerically evaluated (by the system of equations (33)-(34)), the spectrum of generalized 
fractal dimensions Dg (36) can be obtained by the expression for the mass index r(g) (35): 

Np 

In ^ 

D, = ^^ , ■ (38) 

q-1 (g-l)lneL 

Finally, joint using of the Legendre transformation 

a = ^, (39) 
aq 

dr 

/(a) = - ^' (40) 

which sets direct algorithm for transition from the variables {q, r(g)} to the variables {a, /(«)}, and the 
approximate analytical expression (38) for the function Dg makes it possible to determine an expression 
for the multifractal function / (a). 

Now we will consider the special case of search of the box counting dimension Dq and the entropy 
dimension Di. One of goals of this consideration is the validation of assumption of smallness of the 
magnitude D'^^q in the expression (10), which was used for the derivation of log normal distribution of 
the random magnitude pi (13). 

It is easy to show that combined using of Eq. (9) and the inverse Legendre transformation sets, which 
sets transition from the variables {a, f{a)} to the variables {q, 'r(g)}, gives the following dependence of 
a{q) on the moments q: 

a{q) = 2q (ao ~ Do) + ao. (41) 

Substituting (41) and (9) into (37), we obtain the approximate expression for the spectrum of gener- 
alized fractal dimensions Dg(q = 0.1) depending on ao and Dq: 

Dg = ^ [q^ (ao - Do) + qao - Dq] . (42) 

Thus, we can write down the expressions for the box counting dimension Do and the entropy dimen- 
sion Di depending on and Dq 

Do = Do, (43) 

D, = lim ^ ■ ^(^^ - f (-^^^^ = 2Do - ao. (44) 

Here it is necessary to make a few remarks. It will be recalled that f(ag=i) = Di is the value of fractal 
dimension of that subset of the region £, which makes a most contribution to the statistical sum (36) at 
g = 1. However, by virtue of normalizing condition the statistical sum (36) is equal to unity at g = 1 
and does not depend on the cell size e, on which the region £ is divided. Thus, this most contribution 
also is of order unity. Therefore, in this case (and only in this case!) the probabilities of cell occupation 
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Figure 8. Geometrical method for determination of the entropy dimension Di, which leads 
to the obvious equality Di = 2Dq — a^. 




Pi ~ el (6) are inversely proportional to the total number of cells n{e) ~ C^)' the condition 

/(a) = a is fulfilled. 

So, the parameters of system of the equations (33)-(34) obtained by the expression (9) can not in 
essence contain information about the generalized fractal dimensions Dq for absolute value of the mo- 
ments q greater than unity (i.e., |g| > 1). 

Secondly, it is easy to show that the expression (39) for the entropy dimension Di does not depend 
on concrete type of the function /(a), but is determined by its properties, for example, by symmetry 
/(a) = a, f'{a) = 1 and convexity f"{a) > 0. The geometrical method for determination of the 
entropy dimension Di shown in Figure 8 is simultaneously the geometrical proof of assertion (44). 

Thirdly, the expressions for the entropy dimension Di obtained by parabolic approximation of the 
function /(a) and geometrical method (Figure 8) are equivalent. This means that the magnitude D'^^q 
in the expression (10) is equal to zero. Thus, our assumption of smallness of the magnitude D'^^q in the 
expression (10) is mathematically valid. 

In the end, we note that the knowledge of generalized fractal dimensions Dg, the correlation di- 
mension D2 and especially Di, which describes an information loss rate during multifractal dynamic 
evolution, plays the key role for an understanding of the mechanism of secondary aerosol formation, 
since makes it possible to simulate a scaling structure of an atmospheric PM with well-defined typical 
scales. Returning to the initial problem of distribution of points over the fractal set £, it is possible to say 
that the magnitude Di gives an information necessary for determination of point location in some cell, 
while the correlation dimension D2 determines the probability what a distance between the two randomly 
chosen points is less than e^. In other words, when the relative cell size tends to zero (el — )■ 0), these 
magnitudes are anticorrelated, i.e. the entropy Di decreases, while the multifractal correlation function 
D2 increases. 
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5. Conclusions 

Comparative analysis of different pairs of experimental normalized concentration values of atmo- 
spheric PM elements measured in different regions of the Earth shows a stable linear (on a logarithmic 
scale) correlation (r = 1) dependence on different time scales (from average daily to annual). That 
points to a power law increase of every atmospheric PM element mass (volume) and simultaneously to 
the cause of this increase - the fractal nature of the genesis of atmospheric PM. 

Within the framework of multifractal geometry it is shown that the mass (volume) distribution of the 
atmospheric PM elemental components is a log normal distribution, which on the logarithmic scale with 
respect to the random variable (elemental component mass) is identical to the normal distribution. This 
means that the parameters of the two-dimensional normal distribution with respect to corresponding at- 
mospheric PM-multifractal elemental components, which are measured in different regions, are a priory 
connected by equations of direct and inverse linear regression, and the experimental manifestation of this 
fact is the linear (on a logarithmic scale) correlation between the concentrations of the same elemental 
components in different sets of experimental atmospheric PM data. 

We would like to note here that a degree of our understanding of the mechanism of atmospheric PM 
formation, which due to aggregation on inoculating centres (primary aerosols (Dp < l/xm)) show a 
scaling structure with well-defined typical scales, can be described by the known phrase: "... we do not 
know till now why clusters become fractals, however we begin to understand how their fractal structure 
is realized, and how their fractal dimension is related to the physical process" [23]. This made it possible 
to show that the spectrum of fractal dimensions of multifractal, which is a multicomponent (by elements) 
aerosol, always predetermines the log normal type of statistics or, in other words, the log normal type of 
mass (volume) distribution of the i-th component of atmospheric PM. 

It is theoretically shown, how solving the system of nonlinear equations composed of the first mo- 
ments (the average and variance) of a log normal and normal distributions, it is possible to determine 
the multifractal function f(a) and spectrum of fractal dimensions Dg for separate averaged atmospheric 
PM, which are the global characteristics of genesis of atmospheric PM and does not depend on the local 
place of registration (measurement). 

We should note here that the results of this work allow an approach to formulation of the very impor- 
tant problem of aerosol dynamics and its implications for global aerosol climatology, which is connected 
with the global atmospheric circulation and the life cycle of troposphere aerosols [3,21]. It is known that 
absorption by the Earth's solar short-wave radiation at the given point is not compensated by outgoing 
long-wave radiation, although the integral heat balance is constant. This constant is supported by trans- 
fer of excess tropical heat energy to high-latitude regions by the aid of natural oceanic and atmospheric 
transport, which provides the stable heat regime of the Earth. It is evident that using data about elemental 
and dispersed atmospheric PM composition in different regions of the Earth which are "broader-based" 
than today, one can create the map of latitude atmospheric PM mass and size distribution. This would 
allow an analysis of the interconnection between processes of ocean-atmosphere circulation and atmo- 
spheric PM genesis through the surprising ability of atmospheric PM for long range transfer, in spite of 
its short "lifetime" (about 10 days) in the troposphere. If also to take into consideration the evident pos- 
sibility of determination of latitude inoculating centers (i.e., primary aerosol) distribution, this can lead 
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to a deeper understanding of the details of aerosol formation and evolution, since the natural heat and dy- 
namic oscillations of the global ocean and atmosphere are quite significant and should impact influence 
primary aerosol formation dynamics and fractal genesis of secondary atmospheric aerosol, respectively. 

It is important to note also that continuous monitoring of the main characteristics of South Pole 
aerosols as a standard of relatively pure air, and the aerosols of large cities, which are powerful sources 
of anthropogenic pollution, allows to determine the change of chemical and dispersed compositions of 
aerosol pollution. Such data are necessary for a scientifically-founded health evaluation of environmental 
quality, as well as for the planning and development of an air pollution decrease strategy in cities. 
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